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ABSTRACT 


In this paper, a deterministic model incorporating climate changes is formulated 
and analysed to study the impact of temperature factor and variation of river 
water levels on the Schistosomiasis disease transmission in human. Qualitative 
analysis is carried out to determine the basic reproduction number and 
equilibrium points. It is found that the disease free equilibrium is locally 
asymptotically stable if the reproduction number is less unity and endemic 
equilibrium point is globally asymptotically stable if the reproduction number is 
greater than unity. Numerical simulations are performed and suggest that at a 
certain range of temperature, i.e. 15.6° C to 30° C, the rate of humans and snails 
infection of the disease is high. On the other hand, the results show that the force 
of infection experienced by river water levels is high when the river water levels 
is low. The study therefore, recommends that treatments, good sanitation 
practices and human behavioural changes which reduce contact with water 
bodies should be practiced and applied in the community so as to reduce the 
infection intensity in the population. 


Keywords: Climate Changes, Schistosomiasis, Epidemic in Human, Impact, 
Temperature 


INTRODUCTION 


Schistosomiasis, also known as snail fever or bilharzia, is a disease caused by 
parasitic flatworms called schistosomes (Colley et al., 2014). The urinary tract or 
the intestines may be infected. Symptoms include abdominal pain, diarrhoea, 
bloody stool, or blood in the urine. People are infected during routine 
agricultural, domestic, occupational, and recreational activities, which expose 
them to infested water. Lack of hygiene and certain play habits of school-aged 
children such as swimming or fishing in infested water make them vulnerable to 
infection. 

In recent years, the impact of climate changes on human health has attracted 
considerable attention to many researchers. The effects on schistosomiasis 
epidemic diseases have been of an exceptional interest because of its disease 
burden and its transmission sensitivity to environmental conditions. The 
atmosphere, biosphere and hydrosphere are being perturbed in a manner which 
may prove irreversible (Martens et al., 1995). Different studies have shown that 
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the transmission potential of vector-borne diseases is very sensitive to climate changes. Schistosomiasis is the human infection 
probably with the most complex biological cycle, which involves at least two host species (human as a definitive host and snail as 
an intermediate host), two free living transmission stages of the parasite (cercariae and miracidia) and distinct environment (Yang, 
2003). 

Hu et al., (2015), investigated the spatio-temporal variations of S. japonicum infection risk in Anhui province and assessed the 
association of the disease with key environmental factors with the aim of understanding the mechanism of the disease and seeking 
clues to effective control of the disease. Their results suggested that seasonal variation of the normalized difference vegetation 
index, variation of land surface temperature at daytime and distance to the river were negatively significantly associated with the 
risk of schistosomiasis. So far not much has been done on how temperature factor and volume of river water affect the transmission 
dynamics of the disease in human. Therefore, in this paper, it is intended to formulate and analyse a mathematical model that 
determines the effects of temperature factors and variation of river water levels on the schistosomiasis disease transmission in 


human. 


2. MODEL FORMULATION 


In this section, a basic deterministic schistosomiasis model which explicitly incorporates environmental factors such as the effect of 
river water level where miracidia, cercariae and snails live and the effect of temperature in human is formulated and analysed. The 
model involves humans and snails populations respectively. 

The model that is to be formulated is intended to analyse the impact of climate changes on the level of river water where cercariae, 


miracidia and snails live and the temperature factor on schistosomiasis transmission to humans. The model consists of six 


opulations at any time be namely: Susceptible humans S, (ft , infected humans Li , treated humans T(t , recovered 
pop y y P 1 1 


humans R (t ) , susceptible snails S 2 (t ) and infected snails [ 2 (t ) : 

In this model, at any time t , hew susceptible humans and snails recruits enter the populations through birth at constant rates of 

4, and 7, respectively. There is a constant human natural death rate {, and snails natural death rate {,. Infected humans and 

snails die due to disease at constant rates of ™, and ™, respectively. Similarly, infected humans recover from the disease or 

infection due to treatment at a rate of 7, as well as natural recovery at arate of y. N, (t ) is the total humans population given by 
N,(t)=S,(t)+1,(t)+T(t)+R(t). (1) 

Susceptible humans acquire schistosomiasis through infection by infected snails in the physical water environment at a rate of 


A, (t) where 


a(y- Be 2) 


with B, (T) being the rate of disease transmission to the susceptible humans by the infected snails in a river which is favoured by 


suitable environmental temperature T; € is the rate of disease transmission by infected snails due to water level and q@ is an 


impact of water level on the disease transmission. 


On the other hand, NV. 2 (t ) is the total intermediate host snails population given by 


N,(t)=S,(t)+ L(t). (3) 
Similarly, susceptible snails acquire schistosomiasis through infection by infected humans in the physical water environment at a 


rate of Ag (t ) where 


Bp (t ) being the rate of disease transmission to the susceptible snails by the infected human waste (urine or feces) in rivers which 


are favoured by suitable environmental temperature T', while b being the rate of disease transmission by infected human due to 
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water level. The description of model state variables and parameters which are to be used in the model formulation are given in the 
table below. 


Table 1: Description of State Variables and Parameters Used in the Model Formulation 


Variable / Parameter Description 
Ss, (t ) The number of susceptible humans at time f 
fi (t ) The number of infected humans at time f¢ 
T (t ) The number of treated individuals at time ¢ 
R (t ) The number of recovered individuals at time f 
55 (t ) The number of susceptible snails at time ¢ 
I, (t ) The number of infected snails at time ¢ 
Y Natural recovery rate of infected humans 
B (T) Disease transmission rate to the susceptible humans by the 
: infected snails in the river 
B (T) Disease transmission rate to the susceptible snails by the 
. infected human waste in the river 
r, Recovery rate of infected humans due to treatment 
qr, Recruitment rate of susceptible humans 
1, Recruitment rate of susceptible snails 
L, Human natural death rate 
LH, Snail natural death rate 
m, Human death rate due to disease 
m, Snail death rate due to disease 
Rate of disease transmission by infected snails due to water 
E 
level 
; Rate of disease transmission by infected humans due to 
water level 
a Impact of water level on disease transmission 


In formulating the model, the following assumptions are taken into consideration: 

1. The transmission of the disease in the human population is only through contact with infected snails in the physical water 
environment; Likewise the transmission of the disease in the snails population is only through contact with infected human 
waste in the physical water environment, 

2. There is no vertical transmission of the disease, 

3. There is no immigration of infectious humans in the population, 

4. There is no immunity in both humans and snails populations, 

5. Infected humans can recover naturally from the disease, 

6. Infected humans recover from the infection due to treatment, 


7. Variations of temperature due to climate changes affect snails population and contact patterns. 


Taking into account the above considerations and assumptions, we have the following compartmental flow diagram: 
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Figure 1: Compartmental Representation of the Model 


From the above flow diagram, the dynamics of the disease is governed by the following system of nonlinear ordinary differential 


equations: 
Ft =m - AS, - 145 
AAAs, FL 7h —(m +4) h 
<= rl, — wT —1T 
Sahar " 
< SA 85 ples 


dl 
ao —(m, + M,)1,, 


where A, is a force of infection of the disease from the snails to humans, given by 


: caer reer (6) 
and A, is also a force of infection of the disease from humans waste to snails, given by 


A, meee cans (7) 
aS, +bl, 


3. MODEL ANALYSIS 


In this section, the model system (5) will be qualitatively analysed to get insight to its disease transmission dynamics and 
understand the effect of climate changes on Schistosomiasis Epidemic in in humans. Threshold parameter which governs 
elimination or persistence of the disease in humans will be determined and analysed. 


3.1. Positivity of Solutions 


We prove that all solutions of the model system (5) with positive initial conditions remain positive for all the time f > 0. This is 
established by the following Lemma: 
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Lemma 1 Let the initial data be {(5,(0).1,(0).7 (0), R(0),S, (0),1,(0)) 20h eQ. Then, the solution set 


{S,(¢).1,(¢),T (2), R(t),S>(t),Z(2)} of the model system (5) is positive for all f 2 0. 


Proof: From the first equation of the model system (5), we have 


dS 
rr =7,-A8,— 4,5, 
or 
dS 
Fie +u)s, @ 


The inequality (8) has a solution 
f(a +14 Jat 
S,(t)=S,(O)e ° a0) 


where 


ALD 


aS,+el,- 


Similarly, it can be shown that the remaining equations of the model system (5) have positive solutions. 
Hence as f —> ©, the solution set {S, (2),J, (t), T(t), R(t), De (t),J, (r)} of the model system (5) is positive for all f 2 0. 


3.2. Invariant Region 
Since the model (5) is Schistosomiasis model dealing with humans and snails populations, we assume that all state variables and 


parameters of the model are positive for all f 2 0. The model (5) will then be analysed in suitable feasible region where all state 


variables are uniformly bounded in a proper subset QAcII ; 


The total population sizes of human N , and snails N. 4 are given by 
N,=S,+1,+T+R 

and 
N,=S,+1,. 


respectively. 


From the differential equations of the total human population, we have 


dN, dS, dl, dT dR 
— = — + — + — + —, 
dt dt dt dt dt 
=7,-(S,+1,+T+R)u4,-ml,, 


=1,-Nu,-ml,. (9) 


It is noted that in the absence of schistosomiasis, there are no deaths of humans due to the disease, that is 7, = 0. Consequently, 


equation (9) becomes 


dN, 
—<z7,-Nu. 10 
A 1 4 (10) 


The inequality (10) has a solution of the form 
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= —N, (0) 4 


N,(t) ett (11) 
Fadl vad 
nae a 
As t —>00, the population size N, (t ) approaches — . 
Mh 
Similarly, for the population of snails we have 
dN, dS, dl 
Ser (12) 
dt dt dt 
=H, =(S; +E) —mt,, 
= f= Naf, =i T 5. (13) 


In the absence of schistosomiasis disease, there are no deaths of the snails from the disease, that is m, = 0, implying that 
dN, 
dt 


By = Nf (14) 


The inequality (14) has a solution of the form 
Hy | — N, Oe oll 


N,(t)< 
a(t) 120) roo) 


(15) 


pee uD 
As t —>00, the population size NV, (t ) approaches —. 
Hy 


It can be seen that, the humans and snails populations tend to zero as f —> 00. It is therefore concluded that, all feasible solutions 
of the model system (5) are positive and eventually enter the invariant attracting region 


Q=4(S,,1,,7,R,S,1,)€0° 0<N,<,0<N, <=}. (16) 
fal lao) 


Thus, whenever 77, = ft, and 7, = by, Q) is positively invariant and attracting. Then the basic mathematical model is well posed 


and it is mathematically relevant, therefore it is sufficient to study the dynamics of the epidemiological system (5) in the region Q.. 


3.3. Existence and Stability of Equilibrium Points 
In this section, the existence of equilibrium points of the model system (5) and their stabilities are determined. 


Let E (s : we i ; ‘i ,R : JS 7 wi i ) be the equilibrium points of the model system (5). Then the equilibrium points of the model system 


are obtained by setting the right hand side of the model system (5) equal to zero, that is: 
O=7,-A,S; — 4,5; 
O=A4S) —1I, — 71, —(m + 4) 1 
O=nl, — T° -7T 
O=yl,+rT -u4R (17) 
0=7,-A,S; - 4,8, 
0=A,S,—(m, + Wy) 1; 
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Note that, the population of both humans and snails will never be extinct as long as humans and snails recruitment terms 7, and 


7, are not zero respectively. This implies that there is no trivial equilibrium point. 


3.4. Disease Free-Equilibrium Point, E,, 
Disease free equilibrium point (DFE), is a steady state solution of the system when there is no disease. Therefore, in the absence of 


infection, that is, when [ i =T[ = 0), the model system (17) will reduce to 


O=7,- 45, 

which implies that 
* a 
S,=—. (18) 
My 

Similarly 

O=-(u4,+n)T, 
that is, 

T’ =0. (19) 
Also 

O=rT- wR. 
But T= 0 ; therefore R'=0. 
Furthermore, 

O=2,— 15S), 
that is 

* 1 
s,=—. (20) 
Hy 

Thus, from the equations (17), the disease free equilibrium point of the model is given by 

E, =(S/,0, T°,R",S;,0)=|—, 0, 0, 0, —, 0}. (21) 


fa Hy 


3.5. The Basic Reproduction Number, R,, 


Diekmann et al., (1990), defined the basic reproduction number K,, as the average number of secondary infections produced by a 


single infectious host, during his or her entire period of infectiousness. The basic reproduction number is one of the most important 
non-dimensional tool in epidemiology as it sets the threshold in the study of a disease both for predicting its outbreak and for 


evaluating its control strategies. 


The basic reproduction number R,, is computed by using the next generation operator method (Van den Driessche and 


Watmough, 2002) as follows: 
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Let 
(i) F, (x ) be the rate of appearance of new infection in compartment 1, 


ii) V.” (x) be the transfer of individuals into compartment i by all other means, 
i P y 


(iii) V,- (x ) be the rate of transfer of individuals out of compartment i where it is assumed that each of the functions F,, V,” and 


L 


V, are continuously differentiable at least twice with respect to each variable and V, =V, — Ve, 


t l 1 


(iv) X 9 be the disease free equilibrium point. 


According to Diekmann and Heestereek (2000), we call FV the next generation matrix for the model and set the reproduction 


number equal to 


Ry =p(FV ) 
where 
-1 
F= a (%) and V= ad (%) 
x; OX, 


with 121, for the number of compartments and 1< j < m for the infected compartments only. p(A) denotes the spectral 

radius of a matrix A. F and V are mxm matrices, where 71 is the number of infected classes. Consider an infected individual 
» 2 \th 

introduced into compartment k ofa disease-free population. The (i, ji) entry of F is the rate at which an infected individual in 


. ; : . . th 1-5 : 2 
compartment J produces new infections in compartment 1, and the ( ik) entry of V_ is the average time an infected 


7 og Nt zs 
individual spends in compartments J during its lifetime in compartment k . Hence, the (i, k) ‘ entry of the product FV ‘is 
the expected number of new infections in compartment 1 produced by the infected individual originally introduced into 


compartment k (Van den Driessche and Watmough, 2002). We re-write the system (5) as 


A= AS; -r1,-71, —(m, +u,)1,, 

“ ass (m, +4, )I,, 

q = Hla HT, et 
Sa rh +n —R, 

Stan, 50) =f Sy 

2 =~ 2S, — jy 


Following Van den Driessche and Watmough (2002), we define KF and V, from the system (22) as 
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Bh, (T)S\, 


aS,+ € 1, 
F = 
B,(T)S1, 


aS, + bI, 


and 
- nl t+yl, +(m, + 1) 1, 
(m, +)1, 


i 


Now the partial derivatives of (23) with respect to [, and I, gives 


; ap,(T)S? | 
wine (aS,+e1,) 
| __aBS3 é 
(aS, +bI,) 


At disease-free equilibrium point E, 5 


> Alt) 
‘|AIT) 


Similarly, the partial differentiation of (24) yields 


py | oemeaen 0 
‘ 0 (m, + 4) 


which when evaluated a E,, gives 


l(i+mtuty) 0 
ovs)-| 0 (m, +4) | 


1 


Note that V = DV, (a) . Now when we take the inverse of matrix (27) we get 


1 


——————— 0 
y= m+h+h+y 
1 
0 
mM, + My 
Consequently, 
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(25) 


(26) 


(27) 


(28) 
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0 


FV'= (29) 
Py (T) 


(y+m, +7, + 1) 


j . é , . 5 ‘ A -l. 
The basic reproduction number, R, is defined as the spectral radius (dominant eigenvalue) of the matrix FV ie. 


Ry =p (FV") is determined by calculating the eigenvalues of the matrix (29) above. The calculated eigenvalues of matrix (29) 


above are given by 


B(T) A(T) 


al(y +4 +7 +m, )(m, + 4) 


A, =- 


It is clearly seen that Ay is the dominant eigenvalue, hence it is the reproduction number of the model. Therefore, 


Ry = _vAMT)AT) (30) 


al(y +m, +7 +m, )(m, + 1, ) 
3.6. Local Stability of the Disease Free-Equilibrium Point 


is . sess : Xr, Ty ' . : : 
The local stability of the disease-free equilibrium point E, =|—, 0, 0, 0, —, 0] is examined by the linearized form of the 
1 HM, 


model system (5) at the given steady state Ea. We linearize the model system (5) by computing its Jacobian matrix J ge We re- 


define the system (5) as 


fF (Spit. &,S,.4,)=4,-48,- 1S, 
AS, - 71-71, —(m, + ,)1, 


)= 
f,(S,1.T.R,S,1,) = 
)=nl,—w Tt —4T, 
)= 
= 
)= 


( 

( 
FAS iil Als 
f,(S,,.0,,T,R,S5,1,)=nl, +47 — wR, (31) 
i\Sulsl Roots 

( 


FAST RSL 


Wy —AgSy — fyS>, 
AgS, — (fy +m, )I,. 


Hence, at the steady states, the Jacobian matrix of f,, f,, 7, »f,f; and f, with respect to S,,/,,7,R,S, and I, is given by 


ISSN 2278-5469 EISSN 2278-5450 | OPEN ACCESS 


Page285 


DISCOVERY | ANALYSIS ARTICLE 


| of, of, of, of, of, of, | 
“(f) “(e) =(£) =~(z£) ~(z£,) =(E 
25, | 0) ar | 0) ar | 0) aR | 0) a5, | 0) aI, | 0) 
of, ap of, of, of, of, 
=i (fF ). —2(8 ) —2(2) 22(8) =2(2) =@ 
25, | 0) a 0) ar | 0) aR 0) a5, | 0) a, | 0) 
of, Of, Of, of, of, of, 
—i(p) 22(F.) =2(8.) =2(2,) 2h) - 2k 
_ 25, | 0) ar | 0) ar | 0) aR 0) a5, | 0) aI, | 0) 
| Of, of, of, of, of, of, 
ate) by (py) en) Sry) ae 
25, | 0) ar | 0) ar | 0) ri ) a5, | 0) aI, | 0) 
Of; of, Of: of, Of; Of; 
itp) 22h) 228 \ 22(8 |) (2) —202 
25, | 0) a 0) ar | 0) aR | 0) a5, | 0) aI, | 0) 
of, Ff, Of of Of, of 
bf) 28(\ 25(p) 2o(R) 28(p) 28(R 
Las, | 0) ar | 0) ar | 0) aR | 0) a5. | 0) ar | 0) 
which turns out to be 
=(F Fe) 0 0 0 0 Fi | 
F,, -F,, 0 0 0 F,, 
j.= 0 i, -(4,+7) 0 0 0 ee 
“| 0 y A Hh 0 0 = 
0 Fi 0 0 =(7 +15) 0 
| 0 Fy 0 0 Fos —(m, + Ly) | 
where 
Fi, = ate +2), Fie = nie a) Fy, = ies =), Fy =[,+m+4,4+y7], 
(aS,+€/,) (aS,+€1,) (aS, +1,) 
rs Bas; F..-|- B,aS; F.=|- B,bly a B,aS; 
* \(aSteL)y | ~ | (aS,+bL,y |” | (a8,+b4)y | ° | (aS, +01, 
and F.. = _ Bybly : 
» (aS, +b1,) 


Now, when we substitute the values of Ey in the matrix (32) above, at disease-free equilibrium point the matrix becomes: 


-u, 0 0 0 0 A. 
a 
U- <5 0 0 0 B. 
a 
0 +r -(mt+nr) 0 0 0 
i= 1 ie (33) 
a a K “4, =O 0 
B 
a 0 0 -n“, 0 
0 - 0 0 0 ~(m, +s) 
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where 


Ly, =[(,+m+m+7)]. 


We compute the eigenvalues A of the matrix (33) above from the determinant J ky | = 0 where 


=a. 0 0 0 0 B, 
a 
6 ina) 0 0 A, 
a 
0 i, —L,,—A 0 0 0 
_ "= a 
i 0 ig i “Mm-A 90 0 
0 0 0 w-Aa oO 
a 
0 Ps 0 0 (a a 
L a 4 
with 
Ly =((7, +m, + Ly +7), L,, = (4, +7), and 66 =(m, =U): 
From the matrix (34) above, it is clearly seen that there are three eigenvalues which have negative real parts i.e., A, =-LL, 
Ay =—,, and A; = —,. The remaining three eigenvalues are obtained from the determinant 
-(K+m+u,t+7)-A 0 By 
a 
i, —(m,+1,)-A 0 =0. (35) 
B 
= 0 — (m, + Lb ) —A 
a 
This gives 
BPs 
(A+y+m,+7,+u,)(A+m, +) 2 =(), (36) 
If we let 
A =(y+m ++), 
A, = (m, + ) > 
and 
C = PP . 
a 
We can write (36) as 
(A, +4)(A,+4)-C =0. (37) 
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Thus 
A’ +(A,+A,)A+A,A,-—C=0. (38) 


According to Routh-Hurwitz criteria, E, is locally asymptotically stable if A, + A, >O and A,A, —C > 0. 


Then it can be shown that 


A, +A, =(m, +n +,+7)+(m, + u,)>0, (39) 
and 
2 _ BP; 
AA, C=(m+7A+u,+7)(m+4) © >0. (40) 
Thus 
PL <(m +5 +m +7)(m +H), 
or 
, LB, <1. (41) 
a’ (m+n +u,+7)(m, +4) 
This yields 
BB. m 
a (m+7+u,+7)(m, + 4) 
implying that 
T T 
R, = B(T)B.(T) <1. (42) 


(m, +1 + Ly +7)(m, + ft, ) 
Thus, E, is locally asymptotically stable if R, <1. 


3.7. The Endemic Equilibrium State and its Stability 
Endemic equilibrium point E, , is a steady state solution when the disease persist in the population. At the endemic equilibrium 


state, both humans and snails in the population are infected by schistosomiasis, and note that the endemic equilibrium is given by 
Ee (S et eS 2, |. (43) 


By considering the model system (5), the endemic equilibrium value for susceptible humans is given by 


* Tl, 


5) = (44) 
Mt” 
The endemic equilibrium value for infected humans with schistosomiasis is given by 
* TT . 
L=s A (45) 
(Ar +44, )(m + 4,474 +7) 
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Likewise the endemic equilibrium point for individuals getting treatment is given by 


Tr. mrA, 


ee (46) 
(A; + 14,)(m, +, +7, +7) (4, +7) 
In the case of recovered class, the endemic equilibrium point is given by 
* 2 
aA (m+, +77) 
R =, TOO _. (47) 
(Al + 14, )(m, + oH +7) (My +) oh 
In case of snails, the endemic equilibrium value for susceptible snails in the population is given by 
* TU 
—— (48) 
My +A, 
Likewise the endemic equilibrium value for infected snails is given by 
. 2S 1A, 
fo pide) (49) 


. M, + My 7 (1 +4; )(m, +) 


3.7.1. The Existence of the Endemic Equilibrium Point, E, 


In the presence of schistosomiasis disease, that is J, #0 and I, #0, the model system (5), has an endemic equilibrium point 
given by 
Be(Sichet oh ik, )e0) (50) 


where 


(S),0,.T R',S,,5,)>0 and S| = 1 =T'=R' =S; =1, =0. 


Let E = (Silt Ro, be a constant solution of the model system (5). We then express E = (S10 ok 5) 


* 
in terms of I, as 


s'(1)= m,(a+e€1,) 
uy (at+el,)+ Bl, 
I (75) = mB, 
(7, + 1,)(y +m, +7 +4 )(LB, +(at+ €1,) 4) 


r (f)= mi Bd 
(+ 4,)(y+m +74 4)(DB +(ateL) Hm) 

R'(In)= 1,8, (7 +77, + 47) 1, | es 
(4+ 14)(y +m, +7, + 4)(1,8,+(a+e1,)“) 

s*( = m,(a+DbA) 
B,A+t (a+bA) 
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i * Al, 
I, )}=——— 
A(t) at+el, 


an a BP, 
Aa eran earim 


where 


A=(x+m)(y+m, +7 +4,)(LB +(atel)u). 
Now we substitute the expression (51) for 1, which is given by 


dl 
ap Sa (th + Ha) (52) 


at the endemic equilibrium we get 


1] A (mt) |=0 (53) 
1B Bil, +ap,A+ L,b1,7,f, 


where 


A=(K+4,)(y+m, +4 +4,)(LB, +(a+el,)u). 


Equation (53) gives [ i =0 which corresponds to the disease free-equilibrium point and at the endemic equilibrium point it can be 


shown that i, >1 by considering the same expression (53) with some assumptions. By considering the equation 


7. 
Al 1 BB, -(m +44) |=0 
BB, + ap,A+ L,b1,7,f, 
at disease free-equilibrium point where I, =O we obtain 
TT. 
; 1 BP, —(m, + 4,)=0 
a (uy +H)(m +H + +7) yee 
BP _ ML, (44, “ i,) > 1 (54) 
a (u,+7)(m+4+4,+7) 1A, 
Thus 
: peace > 1 ; 
a (u, +7, )(m, + My + +7) 
We then have the following expression which belong to the existence of the disease at endemic equilibrium: 
1 
R, = AP, >0. (55) 


a (44 +7 )(m, +, +, +7) 
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We can only deduce from expression (55) that the only one positive endemic equilibrium exists for Ry > 0. Consequently, there 


exist one unique endemic equilibrium for model system (5) whenever R, > 0. 


3.7.2. Global stability of the endemic equilibrium point, E, 


Theorem 4: If R, >1, the endemic equilibrium E, of the model system (5) is globally asymptotically stable otherwise it is 


unstable. 


Proof: 
We develop the global stability of the endemic equilibrium , by forming the following Lyapunov function: 
VS TP 1) defined by 


Ve) SS) 38 ne |e in| Pa 
S, ra T 
{a-a'-r'in }o[s,s-stin& }o[ 4-1-2 Zh) 
R S i 


_— 
I, \dl, dl, {T- T’ \dT 
= T, dt T dt 
Bs dt f dt 
dS, dl, dT dR dS, dl, 
Now, when we substitute the expressions for from the model system (5) to obtain 


dt’ dt dt’ dt’ dt’ dt 


dv _ S,—S, i, BSI, — 1,5, 
dt S, aS, + € 1, 


Jes +77 -,R] (57) 


=k B,Sol, 
+| —— || ———_--(m,+4)I, |}. 
L, | 2 ete) 
Let S, =S,+S,-S,,],=1,-1,,T =T-T' , R=R-R’, S, =S,—S, and I, =1,-I,. Then 
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# (Ses, areas u4(S) s) 


dt S, a(S,-S))+e(1,-F; 
Se Goer pesmi) 
{EF |[a(n-1)-Wen(r-1)] _ 
(EE \(n-n)en(t-1)-m(R-8)] 


(8:8 AMD 6s, 


(S,-S;)+b(7,-1) 


u(t) $C “(ms )(b-8)] 


L, (S,-S;)+b(1,-1)) 


Let P be equal to values of positive terms and Q be equal to values of negative terms of equation (58). Then after the expansion 


and collection of like terms together in the system we obtain 


“=P-9 (59) 
where 
P= (5, - Si) x, ns B,(S, -S')(1, -1)(h, -I') Tals -f;)) 
5, 1,(a(S,-5;)+e(Z,-4)) T 
Rot )en(t-7)] (5-8) os 
R S, 
B,(S,-S;)(h-h \-4) 
+ , 
1,(a(S, -5;)+b(4,-1)} 
and 


gales Alict) | UntYtremeaen 


s, S,-S})+e(L,-f; I 


1 


(T-T") (+7) (R-R'Y 4, 
R (61) 
a B,(1-1) (1,-15)(m, + 4) 
a 


S,  |a(S,-S;)+b(1,-f 
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dV dV d 
Thus, from equations (60) and (61), if P < Q, then rs will be negative definite, meaning that 7 < 0. It follows that ae =) 
t t t 


if and only if Si; = Sy, I, ae T=T , R=R, 83 = SS I,= L, Therefore the largest compact invariant set in 
* ok So oe * co dV “ . “ ; aon “ 
(s; phigh ads 5 S}.15) EQ; ae = (> is the singleton 12} where E, is the endemic equilibrium of the model system (5). 
t 


By LaSalles’s invariance principle, then it implies that E, is globally asymptotically stable in Q if P<Q. 


4. NUMERICAL SIMULATIONS OF THE MODEL 


In this section, we illustrate the analytical results of the study by carrying out numerical simulations of the model system (5) using a 
set of parameter values given in table 2. The values are obtained from published literature and other assumed parameter values 


with the initial conditions: S, (0) = 6500, 1, (0) =40, S,(0) =15000, 7,(0)=300, W(0)=22, T(0) =15.6°C. 


Table 2: Model parameter Values and their interpretations to be used in the Numerical Simulations 


Parameter Meaning Value Source 
1, Recruitment rate of susceptible humans 800 Chiyaka et al., (2009) 
1, Recruitment rate of susceptible snails 2500 Garira et al., (2014) 
Ly Human natural death rate 0.0000384 Feng et al., (2002) 
LH, Snail natural death rate 0.000569 Spear et al.,(2002) 
m, Human death rate due to disease 0.0039 Feng et al.,(2004) 
m, Snail death rate due to disease 0.0004012 Mangal et al., (2008) 
Disease transmission rate to the susceptible 
B, eae 0.406 Spear et al., (2002) 
human by the infected snail in river 
Disease transmission rate to the susceptible 
B ea 0.615 Mangal et al.,(2008) 
snail by the infected human in river 
Z Recovery rate of infected human due to 0.49 Mushayabasa and Bhunu, 
: treatment (2010) 
Q Rate of ground water inflow to the river 15000 UL! Assumed 
Cc Precipitation rate due to rainfall 0.155 Assumed 
n Rate of evaporation 0.66 Assumed 
Rate of environmental increase or decrease 
Oo . . 1.9999 Assumed 
in temperature due to climate change 
Prescribed maximum environmental 
: 99°C Assumed 


ax temperature increase (50—1000°C) 


Rate of disease transmission by infected 
€ : 0.0002 Spear et al.,(2002) 
snail due to water level 


Rate of disease transmission by infected 
0.00022 Assumed 
human due to water level 


Impact of water level on disease 
— 0.00025 Assumed 
transmission 


Figures 2 show the dynamics of state variables for human population. Figure 2(a) shows the susceptible human while figure 2(b) 
shows the infected human in the population. 
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Figures 2: Dynamics of the State Variables of the Schistosomiasis Dynamics on (a) susceptible human and on (b) infected 
human in the population 


From figure 2(a), it is observed that the number of susceptible individuals in the population decreases with time. This happens 
due to susceptible humans who are exposed to infection, die naturally. Similarly, figure 2(b) shows that infected humans increase 
due to higher number of susceptible humans moving from susceptible class to infected class in the population and then decrease 
exponentially time due to schistosomiasis-induced deaths, natural deaths and those who recover from treatment which is practiced 
in the community. 


Figures 3 show the dynamics of susceptible snails and infected snails. 
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Figure 3: State Variables of the Schistosomiasis Dynamics on (a) Susceptible snails and (b) Infected Snails in the Population 


From the two graphs it is seen that the susceptible snails in figure 3(a) are decreasing with time. This is due to the fact that some 
of them are dying naturally while others are moving to the infected class after being infected by infested human wastes in the river. 
Similarly, in figure 3(b) it is observed that there is an exponential increase of infected snails up to a certain level and eventually 
started decreasing with time. The exponential increase of infected snails is due to high number of susceptible snails in the 
population being infected, also the exponential decrease of infected snails with time is due to disease-induced deaths and natural 
death of infected snails. 
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Figures 4 show the rate of environmental increase in temperature due to climate changes and the rate of environmental decrease 


in temperature due to climate changes respectively. 
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Figure 4: Results showing the increase in temperature and the decrease in temperature due to climate changes 


From the figure 4(a) it is observed that as the number of susceptible humans decreases exponentially due to infection. The 
number of infected population is found to be increasing at a certain time interval and then start decreasing with increase of time. 
The other reasons for susceptible class to decrease with increase of time may be is due to natural death, and on the other hand the 
infected class of humans is decreasing perhaps is due to disease-induce rate and natural death in the particular population. On the 
other hand, figure 4(b) shows an addition of treated class. It is seen that the graph of treated individuals is increasing after the 
application of medicine such as praziquantel to the infected population. 


Figures 5 show the evolution with time of infected, treated and recovered population. 
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Figures 5: Results Showing the Evolution of Infected, Treated and Recovered Population with time 
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The simulation results shown by figure 5(a) shows that the high amount of cercariae released by infected snails to the physical 
water environment causes the higher rate of infection at an individual level. The number of infected individuals will continue to 
increase until proper measurements are taken to combat the disease. The decreasing graph of infected groups with increase of time 
indicates that some individuals are recovering due to treatments in the community as shown by treatment graph in the same figure 
5(a). Natural deaths and disease induce death rates are among of the reasons which causes the decrease of infected population. On 
the other hand, figure 5(b) shows the recovered population after proper treatments are taken in the particular population of infected 


individuals. Thus treatment has a positive impact to the society. 


Figures 6 show the numerical solutions of infected, treated and recovered individuals in the population. 
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Figure 6(b) 
Figure 6: Results showing the rate of recruitment for both susceptible human and snails at a particular population 


Figures 6 illustrate the graphs of numerical solutions showing propagation of susceptible humans and susceptible snails for 


different values of rate of recruitments. The rate of recruitment for new susceptible humans are 7, = 80, 200, 400, 600 and 


the rate of recruitment for a new susceptible snails are 7, = 250, 500, 700 and 1000 respectively. We note from figure 6(a) 


that the emerging of new susceptible humans increase the intensity of human population at an individual level, also from figure 
6(b), it is observed that the new susceptible snails increase the intensity of snail population in the physical water environment. 
We note that when the rate of recruitment for susceptible snails is high in the physical water environment and when the rate of 


miracidia released by human wastes increases in the physical water environment lead to the higher number of snails infection rates. 
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The higher the number of infection rate of snails leads to a higher number of cercariae in the physical water environment, which 
upon when in contact with humans accelerates the rate of re-infections again and a big population will likely suffer form the 
disease. It is suggested that elimination of snails by the application of molluscicides is the best solution to reduce or eliminate the 
disease. 


Figures 7 show the natural death rate of susceptible humans and susceptible snails in the physical water environment. 
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Figure 7(b) 
Figures 7: Effects of Variations of River Water Levels on Schistosomiasis Epidemic in Snails 


From the figures 7, it is observed that when the rate of natural death is low, the number of susceptible humans and snails are 
higher in the entire population as indicated by 4, = 0.000384 and “, =0.000014 respectively. As the natural death rate 
increase for both populations, the number of susceptible humans is found to be decreasing at a higher rate as indicated by 
H, = 0.384 . Likewise, the number of susceptible snails in the physical water environment are also removed fast as indicated by 
H, = 0.0014 . It is also noted that the death rate of snails reduce the intensity of disease at individual level, so it is suggested that 


an elimination of snails which act as an intermediate host in the disease transmission in all sources of water bodies will help to 
reduce schistosomiasis in the community. 


Figures 8 show different rates of disease transmission for infected individuals in the population and disease induce death rates. 


ISSN 2278-5469 EISSN 2278-5450 | OPEN ACCESS 


Page297 


DISCOVERY | ANALYSIS ARTICLE 


1400 
1200 
1000 f---- 
c 
2 
3 soo}. 
a 
=) 
a 
3 
= 600} 
2 
= 
400 
200 
i : : : : : : : : 
o O85 1 16 2 25 3 35 4 45 5 
Time (years) 
Figure 8(a) 
1200 
et lam 193609 
: j : : —— m,=0.13699 
1000 + aaa i moe , ere ge na Be ai a aaa m,=0.013699 
= : es m,=0.0013699 
& 
a oe, Se en) Se. See. TRE i ee are eee 
2 
= : 
a : 
oO : 
2 : 
= 600 - 
Ss : 
8 : 
2 E 
= 400 : on 
= : 
Qa 3 
o i L i i i i i i 


in} 0.5 1 AS 2 25 S) 3.6 4 45 5 
Time (years) 


Figure 8(b) 


Figures 8: Dynamics of Infected Humans under Different Rates of Disease Transmission, {, and Disease Induce Death Rate 


mM 


From figure 8(a), it is observed that the rate of disease transmission to the susceptible human by cercariae released by infected 
snails in the physical water environment depends on the ambient temperature. When the ambient temperature is below a certain 


level which do not favour the reproduction of snails, the rate of transmission is also found to be low as shown by B, = 0.028 and 
m, =0.038 above, but when the ambient temperature increases up to a certain level where it become conducive for snails 


reproduction, we observe that the rate of disease transmission increases. This happens due to the fact that more snails get recruited 


in the physical water environment and more snails get infected by miracidia which results in a large amount of cercariae parasites 
which again infects human population at a higher rate as illustrated by B, = (0.058 in the figure 8(a) above. 
The high number of human infection rate increases the death rate due to disease for infected individuals in the population as 


shown by figure 8(b). When the infection is high, the death due to disease is also found to be high as indicated by m, = 0.23699 
and m, = 0.13699. As the rate of transmission decreases with time, it is observed that the number of deaths decreases as shown 


by m = 0.0013699 in figure 8(b). It is suggested that human behaviour changes which reduce contact with infected water 


bodies and good sanitation should be practiced in the community. 


Figure 9 shows the dynamics of infected snails under different rates of disease transmission /3,. 
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Figure 9: Different rates of disease transmission for infected snails in the population 


Figure 9 illustrates graphs of numerical solutions showing different rates of disease transmission to the susceptible snails by the 


infected humans waste in the physical water environment. The rate of snails infection is high when the amount of infected human 
wastes is also large in the physical water environment, and this is shown by B; = 3.50 from the given figure above. As the 
amount of infected waste decreases, it is also observed that the rate of snails infection decreases with increase of time, as seen in the 
figure above where 8, =3.30, £8, =2.90 and 8, =0.6150. Again the transmission of the disease from humans to snails 


depends on the ambient temperature of the environment. The values for different rates of disease transmission to the susceptible 


snails by infected human waste in river are as follows: B; = 3.500, B; = 3.300, 2 = 2.900, B, =0.615. 


Figure 10 shows different values of disease transmission by infected humans due to river water levels b. 


Figure 10 shows the graphs of numerical solutions showing different values of disease transmission by infected human due to 


river water levels. When the levels of river water is low, the rate of disease transmission become high if that amount of water in 
infested by cercariae as it is shown by b=0.0002 in the figure 10. As the levels of river water increases, the force of infection to 
susceptible individuals decreases which leads to the decrease of infection rates for human in the population as shown by b = 0.02 


and b=0.2 respectively. The tested values for rates of disease transmission by infected humans due to infested river water levels 


by cercariae are: b= 0.0002, b=0.002, b=0.02, b=0.2. 


Figure 11 shows different values for disease transmission by infected snails due to water levels, €. 
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Figure 10: Different values of disease transmission by infected human due to water levels 
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Figure 11: different values for disease transmission by infected snails due to water levels 
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From the figure 11, it is observed that number of infected snails which is presented by €= 0.0002 is high only when the levels 
of infected river water by miracidia is low, say €= 0.0002, as the levels of infected river water increases say from €= 0.0002, 


to €=0.002, the number of infected snails are also found to be decreasing with time as seen when €= 0.002. This trend is 
expected to continue until there is no further infection when the levels of river water is maximum. Values for different variations of 


river water levels used in the simulation are: €= 0.0002, e= 0.002, €= 0.02 and €= 0.2. 


Figure 12 shows the recovery rate of infected humans due to different rates of treatments /, . 
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Figure 12: Recovery rates of infected human due to different rates of treatments 


From the figure 12, it is observed that the application of treatments to the infected individuals in the population is an important 
factor in the elimination of the disease. Figure 12, presents recovery rates of infected population by the application of medicine. 


7, =0.82 indicates low amount of medicine supplied to treat the infected people in the community, but as you increase the rate of 
treatments in the population, all infected individuals will start to recover at different rates depending on the rates of treatments 
applied as shown above. % = 0.9 indicates higher rate of recovery because of proper and efficient treatments conducted in the 


community. 


Figures 13 show the variation of natural recovery and recovery rate due to treatment by infected population. 


ISSN 2278-5469 EISSN 2278-5450 | OPEN ACCESS 


PagesO1 


DISCOVERY | ANALYSIS ARTICLE 


Recovered population 


f i i 
0 O46 1 15 2 Pas 3 a5 4 45 5 
Time (years) 


Figure 13(a) 


= =0.0005 
= =0.005 
"950.05 
mmm y.5 


population of recovered human 


Tire (years) 


Figure 13(b) 
Figures 13: Results showing the comparison of natural recovery and recovery due to treatments to humans population 


From figure 13, it is seen that the rate of natural recovery of infected individuals depends on the immunity within an individual. 
If the immunity is strong within an infected human, the rate of natural recovery also become fast as indicated by figure 13(b) for 


vy =0.5 and y =0.05, but the lower the immunity leads to slow natural recovery of individuals as demonstrated by y=0.0005. 
On the other hand the application of strong and proper medicine within an infected individuals helps to cure the disease and 


recover at a higher rate as indicated by figure 13(a) for = 0.9. 


5. CONCLUSIONS AND RECOMMENDATIONS 

In this study, we have formulated and analysed a mathematical model of schistosomiasis concerning the impact of climate changes 
on the disease transmission in human. Temperature variation and the different levels of river water were considered as agents of 
climate changes on the disease transmission. The dynamics of human population and snails population have been incorporated 
with environmental factors such as the effect of climatological environment. Both qualitative and numerical analysis of the model 
have been used to determine the basic reproduction number as well as the existence and stabilities of the model equilibrium points. 
Qualitative analysis showed that the model has two equilibria, namely: the disease free equilibrium point and endemic equilibrium 
point. Using next generation method, it was found that the disease free equilibrium is locally asymptotic stable when the 
reproduction number Ry <0, and unstable otherwise. In the case of endemic equilibrium point, the qualitative analysis was done 


by using Centre manifold theorem developed by Castillo-Chaves and Song (2004). It was concluded that near the threshold 
R, =, there exists a stable endemic equilibrium point which is globally asymptotically stable if R, >0. 
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Numerical simulations were performed and it was observed that the transmission of schistosomiasis in both human and snails 
is effective when environmental temperature favours the reproduction of snails. It was also noted that when the temperature was 
below a certain level, the rate of infection was observed to be decreasing as a results of snails been in the hibernation state. Likewise, 
when the environmental temperature is above a certain level, it was noted that the rate of infection was also decreasing, due to the 
fact that all infected snails become unable to release or shed cercariae to the physical water environment. Furthermore, simulation 
results showed that variation of river water levels where miracidia, cercariae and snails live, had an effect on the disease 
transmission. It was observed that when the levels of river water is low the number of infected humans by schistosomiasis in the 
population is larger than when the levels of river water is high. 

Eradication of schistosomiasis remains a challenge in most of the developing countries in the world, hence there is a need to 
strengthen the control strategies at hand as well as putting more efforts to overcome the epidemic especially to those areas close to 
rivers, dams and ponds. Thus, it is recommended that: 

1. Human behavioural changes which reduce contact with infected water bodies should be practiced so as to reduce the infection 
intensity at individual levels, 

2. Good sanitation practices should be applied by the community which reduce contamination of water bodies so as to reduce the 
rate of infection in humans at individual level, 

3. Elimination of snails population using molluscides should be practiced by public health interventions, 

4. More awareness campaigns on the side-effects of schistosomiasis epidemic to humans, should be practiced especially to those 
areas close to waterbodies. 
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